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The problem of neutrino spin rotation in dense matter and in strong electromagnetic fields is 
solved in accordance with the basic principles of quantum mechanics. We obtain a complete system 
of wave functions for a massive Dirac neutrino with an anomalous magnetic moment which are 
the eigenfunctions of the kinetic momentum operator and have the form of nonspreading wave 
packets. These wave functions enable one to consider the states of neutrino with rotating spin as 
pure quantum states and can be used for calculating probabilities of various processes with the 
neutrino in the framework of the Furry picture. 



INTRODUCTION 

Neutrino physics is one of the most rapidly developing areas of high-energy physics. The fundamental experimental 
result was obtained in this field in recent years — neutrino oscillations were discovered , ll] . From the theory of this 
phenomenon (see, for example, f3|), which is based on the ideas of Pontecorvo and Maki, Nakagawa and Sakata [H , 
it follows that oscillations are possible only when the neutrino mass is nonzero. As a consequence of this circumstance, 
the possibility exists of neutrino spin rotation, i.e. of transitions between the active left-handed neutrino polarization 
state and the sterile right-handed state. As it is understood today, in contrast to flavor oscillations, this process can 
never be observed if the particle moves in vacuum. For the effect of spin rotation to exist some external influence leading 
to effective breaking of the Lorentz-symmetry of the theory is needed. For instance, dense matter or electromagnetic 
fields can serve as the physical reason for this phenomenon. In the latter case, nontrivial electromagnetic properties 
of the neutrino, in particular, the anomalous magnetic moment of a Dirac neutrino can lead to rotation of its spin. 

The conventional approach Q to the theory of neutrino spin rotation, based on solving the Cauchy problem 
for the Schrodinger-type equation with an effective Hamiltonian, was developed for the description of ultrarelativistic 
neutrinos (the particles that are observed in experiments now available). This approach is in fact quasiclassical and 
is equivalent to solving the Bargmann-Michel-Telegdi (BMT) equation in the spinor representation (see @). 

The aim of this work is to develop a consistent quantum theory of neutrino spin rotation to describe neutrinos 
with nonzero mass, including the low-energy region that may play a significant role in astrophysics [lo| . Following 
papers , we use the method of effective potential [ll| that allows us to take into account collective influence of 

background particles on the propagation of neutrino. 

To this end, we find a complete system of solutions of the Dirac-type equation that describe the action of both 
electromagnetic field and dense matter on neutrino dynamics. These solutions have the form of nonspreading wave 
packets. They enable one to consider the states of neutrino with rotating spin as pure quantum states and to evaluate 
probabilities of various processes with participation of a neutrino in the framework of the Furry picture . 

The paper is organized as follows. In the first Section we present the model we are working with and briefly describe 
how the effective potentials are constructed. In Section HI] we summarize the main ideas of the method used for solving 
the Dirac equation. In Section IlIII we introduce a kinetic momentum operator and justify the advantages of using its 
eigenvalues as quantum numbers, instead of the eigenvalues of the canonical momentum operator. Then we find both 
stationary and nonstationary wave functions. Here we also derive the dispersion law and calculate the group velocity 
for a neutrino in the dense magnetized matter. We demonstrate that the group velocity of a neutrino is independent 
of its spin orientation. In Section IIVI we find the explicit expressions for the kinetic momentum and spin projection 
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operators. In SectionlVlwe use the obtained results for quasiclassical interpretation of neutrino spin behavior in dense 
magnetized matter. In the conclusive Section [VII we compare our results with those obtained by other authors. 

The units h = c — 1 are used throughout the paper; metric tensor, = diag(l, —1, —1, —1); totally antisymmetric 
Levi-Civita tensor, e'^'^^^ = +1 if {/ii^pA} is an even permutation of {0, 1, 2, 3}, e^'^P^ = —1 if it is an odd permutation, 
g/ii'pA ^ Q otherwise. 



I. WAVE EQUATION FOR NEUTRINO IN DENSE MATTER AND IN ELECTROMAGNETIC FIELD 

Let us first discuss the model we use for the description of the neutrino behavior in the matter and in the electro- 
magnetic field. At present, there are no experimental data that might help to make the choice between the Dirac or 
the Majorana nature of a neutrino [l3l |. Therefore, in this paper we discuss spin rotation of a massive Dirac neutrino 
with anomalous magnetic moment propagating through dense magnetized matter. 

As a first approximation, we can consider spin rotation independently of variations of the particle fiavor, i.e. we can 
deal with the mass eigenstates of neutrino. It is quite obvious, due to the fact that the physical origins of neutrino 
flavor oscillations and its spin rotation are different. The Dirac-Pauli equation [3] serves as the basis of this approach 
in electromagnetic fields. Just in this way, the behavior of the neutrino spin in a constant homogeneous magnetic field 
was investigated in pLS;] . In order to consider the influence of matter on neutrino dynamics through weak interaction, 
the Dirac-Pauli equation should be modified. 

As it was proposed in (TTj . [l6j . if matter density is high enough for considering weak interaction of neutrino 
with the background fermions as coherent, it is possible to describe neutrino interaction with matter by an effective 
potential. To get covariant realization of this idea, it is necessary to take into account that the term that describes 
direct interaction of neutrino with matter must depend on the currents 



jf = {nfU%nfUf}, (1.1) 



and on the polarizations 



of the background fermions / |17H19l | . It is significant that only jj and X'j characterize matter integrally. In formulas 
(jl.ip and (|1.2p . Uf and C/ (0 ^ IC/P ^ 1) a-re the number density and the mean value of the polarization vector of 
the background fermions / in the center-of-mass system of matter (w^ ~ {u'j,Uf} denotes the four-velocity of the 
reference frame), where the mean momentum of the fermions / is equal to zero. 

In the framework of the standard model (for a more general case, see Appendix [X)) . the interaction term can be 
derived directly from the model Lagrangian. To this end, it is necessary to write this Lagrangian in the low-energy 
effective form using the four-fermion interaction and replace the vector ipf(x)j^ipf{x) and axial vector tpf(x)j^^^tpfi^) 
currents of the background fermions by and respectively. The explicit form of jj and A^ is obtained as a 
result of averaging of currents over the fermion statistical distribution function and depends on the type of the matter 
considered. In particular, the expression for the polarization A^ of electron gas in a magnetic field at finite temperature 
can be found in Note that, since the origin of the effective potential is forward elastic scattering of neutrino, it 
can also be calculated using field-theoretic methods po| . 

As a result we come to the equation 



- T^^'^Ua + 7') - T;^^oF'^''^^^^ - rn ] ^{x) - 0, (1.3) 



where i^'^" is the electromagnetic field tensor. The 7-matrices satisfy 7^7'' -f 7''7'^ = Sg^'' with 7° Hermitian, 7 anti- 
Hermitian, and are related to the /3 and a matrices through 7° = /3,7 = /3a; a^^ — 5(7^7"— 7^7^), 7^ = — «7'^7^7^7'^. 
The effective four-potential is a linear combination of currents and polarizations of the background fermions 

r-Hipfff^pf^'f)- (1-4) 
/ 

Summation in (11.41) is carried out over all fermions / of matter. The expressions for the coefficients Py^'^'* are 
determined as 

pf=V2GF[h, + T^/^-2Q^f^sin^e^Y pf ^ -^Gy [^h, + T^^^Y (1.5) 
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Here Q^-''^ is the electric charge of the fermion /; t!^'^' is the third component of the weak isospin; Gp and 6*^ are 
the Fermi constant and the Weinberg angle respectively; I^iy = 1 for the electron neutrino interaction with electrons, 
I^i, = — 1 for interaction with positrons, otherwise = 0. 

While equation (|1.3I) describes mass states of neutrino, potential (|1.4p is flavor dependent, and in the general 
case, this leads to correlations between spin rotation and flavor oscillations. In order to construct a mathematically 
consistent approach to a description of neutrino spin rotation, we have to avoid the correlations. This is possible if we 
assume that the effective potential which describes the influence of dense matter on neutrino is the same for different 
flavors. In this case we can construct flavor states of the neutrino as the linear combination of its mass states with 
coefficients that are elements of the mixing matrix of the neutrino in vacuum. 

The physical justification of this model can be made in the following way. In order to consider effective potentials 
as flavor independent, it is necessary to assume that the fraction of electrons in the matter is small, i.e. the electron 
number density rig is approximately equal to zero. Calculations which are in a satisfactory agreement with the 
experimental data show that in the center of a neutron star the fraction of electrons does not exceed a few per cent 
[2l|. Therefore the model approximation — is quite appropriate. 

We restrict ourselves to considering equation p.3p in the case of a constant homogeneous electromagnetic field and 
constant currents and polarizations of matter, 

F^" = const, = const, = const, (1.6) 

as the first order approximation of the realistic background. Because of p.6p additional conditions for F^'^ and 
may be obtained. The strengths of the electric and magnetic fields and average currents and polarizations of the 
background particles should obey the self-consistent system of equations including the Maxwell equations, the Lorentz 
equation 

and the the BMT quasiclassical spin evolution equation 

X}. (1.8) 

Here the dot denotes differentiation with respect to the proper times r/ of particles. From equations (|1.7p and (|1.8p 
we find that for charged particles the conditions = 0, = are equivalent to F'^^jj = 0, F^^A^ = 0. However, 
it seems reasonable to assume that in dense matter velocities of the center-of-mass systems for all components are 
equal, thus for the neutral particles similar conditions should hold as well. In this way we obtain the restriction 

F^^f^ = 0. (1.9) 

It should be emphasized that condition (|1.9|) is the result of the fact that average currents and polarizations of particles 
of matter in an external field should satisfy the classical equations of motion. The physical meaning of this condition 
is discussed in more detail in Section fVl 

Note that equation (|1.3p with constant coefficients appears in the standard model extension (2^, as well. In 
this case, F'^'^ and describe vacuum condensates that break the Lorentz-symmetry of the theory. 

II. FORMULATION OF THE PROBLEM 

The Dirac equation, in particular equation (|1.3p . is a partial differential equation. Therefore, as it is known, its 
general solution is defined up to an arbitrary function. However, in quantum-mechanical applications we deal with 
the so-called complete integral which depends on a set of constants, i.e. quantum numbers. Those are the eigenvalues 
of some self-adjoint operators. For the classification of particle states it is necessary to introduce the complete set 
of the operators — integrals of motion. Note that in nonrelativistic quantum mechanics any self-adjoint operator 
may serve as an operator of observable, which is not the case in the relativistic mechanics, where only integrals of 
motion, commuting with the operator of the equation, can be treated as operators of observables [23]. The choice of 
the complete set of them is different in each particular case and should be adequate to the problem being solved. 

In the case of the problem of spin evolution, operators of kinetic momentum components should be included in the 
complete set. It becomes obvious if we start with the following argument. The direction of particle polarization is well 
defined in its rest frame, and then one finds the polarization in the laboratory frame upon carrying out an appropriate 
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Lorentz transformation. This transformation is defined by the group vefocity of the particle, in other words, by its 
kinetic momentum. This imphes that it is not the canonical momentum but the kinetic momentum operator that 
defines the direction of particle propagation. 

Now, there is only one problem that remains to be solved. The form of the kinetic momentum operator for a 
particle with spin propagating under the infiuence of external fields is not known beforehand. So we have to find a 
self-adjoint operator with the eigenvalues q'', which satisfies the condition = and may be interpreted as the 
components of the particle kinetic momentum. 

Let us discuss this issue in more detail. In the mathematical apparatus of quantum field theory, a particle is usually 
identified with an irreducible unitary representation of the Poincare group 25j. The irreducible representations are 
characterized by two invariants of the group 



= W^'W^,^ -m^s{s + l). (2.2) 

The translation generators P'' are identified with the components of the particle momentum, and the Pauli-Lubanski- 
Bargmann vector 

= -]^e^'''P^M,pPx, (2.3) 

where M^'^ are the Lorentz generators, characterizes the particle spin. The invariant is the particle mass squared 
and s is the value of its spin. 

A space of unitary representation is defined by the condition called "the wave equation for a particle with mass m 
and spin s." The wave equation for particles with spin s = 1/2 is the Dirac equation 

(z7'^5^-TO)<f'o(a;) =0. (2.4) 

In this case the realization of generators of the Poincare group and the Pauli-Lubanski-Bargmann vector in the 
coordinate representation is 

p^' = i^^', to'"^ = i(x^5'^ - x'^a^) + ^c^^^ u;^ = ^ 7V^''5^. (2.5) 

These operators commute with the operator of the Dirac equation and can be identified with observables. They 
have a self-adjoint extension on the subsets of solutions of equation (|2.4p with a fixed sign of the energy with regard 
to the standard scalar product, 

{<P,<F) = Jdx<p\x,t)>I'{x,t). (2.6) 

Three-dimensional spin projection operator s{p) is a set of coefficients Si{p) of the expansion of the vector in 
spacelike unit vectors S-'{p) {i = 1, 2, 3): 

s.{p) = - l^w^S>^{p), (2.7) 

where 



p^5f(p) = 0, S':{p)SM^-5^3■ (2-8) 

Obviously, 

[s,(p),Sj(p)] = ie,jkSk{p)- (2.9) 

The choice of these unit vectors is not unique, and it is possible to construct operators that determine the spin 
projection on any direction in an arbitrary Lorentz frame. 

The above description of the particle characteristics cannot be directly used in the presence of external fields, where 
operators (|2.5p are not necessarily integrals of motion. In this case the classification of particle states is usually realized 
by linear combinations of operators p^^ and with coefficients depending on coordinates [26l |. Unfortunately, the 
physical meaning of these operators is usually not quite clear. Therefore, one should formulate a description of the 
particle motion in an external field in the same clear and detailed way as for a free particle. 
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This can be achieved basing on the following ideas. Since an irreducible representation of group is defined accurately 
up to an equivalence transformation, it is reasonable to state the problem of finding such realization of the Lie algebra 
of the Poincare group for which the condition of irreducibility of the representation leads to wave equation describing 
a particle in a given external background. To solve this problem it is necessary to find a unitary operator U{x,xq) 
which converts solutions of the wave equation for a free particle \2A\ to solutions ^{x) of (|1.3I) : 

U{x,xo)Mx)=^{x). (2.10) 

Thus, ?7(x,xo) is an intertwining operator in the sense of Darboux 27]. This operator in our case should satisfy the 
equation 

[^rd^ - ^7Vm(1 + 7') - ^A^oi^'^V^, j [/(x, xo) - U{x, Xq) (^7'^^^) = 0. (2.11) 

Therefore, operators 

= [/(x,xo)p''C/"^(x,xo), m^'' = [/(a;,a;o)m^''C/"i(a;,a;o) (2.12) 

commute with the operator of the wave equation. As a consequence, the Pauli-Lubanski-Bargmann vector and 
the components of the three-dimensional spin projection operator &i can be constructed in the same way as in the 
case of a free particle: 

= - i e'^'^f^m.pPA, 6. = - ^!21J,,^f (p). (2.13) 

The above statement may be reduced to the following: the wave function of a neutrino in dense matter can be 
derived with the help of a solution of the Dirac equation for a free particle and of some unitary evolution operator. 
A complete set of integrals of motion may be constructed with the help of operators (|2.12p . The physical meaning of 
eigenvalues of observables, i.e. quantum numbers, is clear enough then. However, in the general case, U{x,Xo) is an 
integral operator and, as a consequence, operators (|2.12p are also integral ones, so it is difficult to find U{x,xo) by 
direct calculations. That is why we are trying to find the evolution operator for equation (jl.Sp using the correspondence 
principle. 



III. WAVE FUNCTIONS 



When motion of a massive particle with spin in electromagnetic fields is described in the framework of the quasi- 
classical approach, we use two four- vectors, namely the four-velocity and the spin vector obeying the conditions 

u^ = \, (US') = u^S*^ = 0. (3.1) 

These vectors are solutions of the Lorentz and BMT equations, respectively. As it was shown in [285, both four- 
vectors can be constructed with the help of one and the same evolution operator that acts on different initial values of 
these four- vectors, satisfying relations (13.11) . In other words, the evolution operator for the BMT equation completely 
describes quasiclassical behavior of the particle. 

As far as the neutrino is a neutral particle, the BMT equation for the spin vector of a neutrino moving with the 
four-velocity has the form 

S^^2^lo{9^''~u^^u'')F^,Sr (3.2) 

We can extend the BMT equation to include weak interaction between the neutrino and dense matter. It was shown 
phenomenologically in [29] that if effects of neutrino weak interactions are taken into account, the Lorentz- invariant 
generalization of the BMT equation is 

= 2 (5^" - u^^u^) ifioF^, + G„,) ^^ (3.3) 

where 

G'^'' - ^e^'-^VpUA. (3.4) 
We can get an analogous result by averaging the equations of motion for the Heisenberg operators (see [1^ ) . 
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A remarkable feature of the BMT equation is its universality — at least in the zeroth order of the Planck constant, 
particles with an arbitrary spin may be described by the same BMT equation. However, a part of the information 
about the particle behavior gets lost due to distinctions in transformation properties of the states of the particles with 
different spins. This information may be restored in the following way. We can describe a particle using quasiclassical 
spin wave functions ^{t) constructed in such a way that the vector current 

f^iT)=§iT)rnr), (3.5) 

as well as the axial current 

j^(r) = ^iT)^^rnr) (3.6) 

built on their base, obey the BMT equation, and, as a consequence, may be interpreted as four-velocity u'^ and spin 
vector S*^ respectively. 

We can introduce quasiclassical spin wave functions as follows [3l|, [13] ■ Let the Lorentz equation be solved, i.e. the 
dependence of the particle coordinates on proper time is found. Then the BMT equation transforms to an ordinary 
differential equation, whose resolvent determines a one-parameter subgroup of the Lorentz group. The quasiclassical 
spin wave function is a spin tensor, whose evolution is determined by the same one-parameter subgroup. It is easy 
to verify (see Appendix [C| that for a neutral particle, represented by a Dirac bispinor, the equation that defines the 
evolution operator for the wave function ^{t) takes the form 

U{t, to) = -j^a^^if^^q" U{t, to). (3.7) 

TO 



Here (ys^ = f /2 + Hf^^q^/m, H'"' = no*F'"', where *F'"' = -^e'^'P^Fpx is the dual electromagnetic field tensor and 
g'' = mu''. 

Let us look for a solution of equation ()1.3p in the form 

^{x) = e-'-^(^)C/ (r(x), 0) <fo(a;), (3.8) 

where U{t, tq) is an evolution operator for the quasiclassical spin wave function, e"*^^^-' is a phase factor, and ^^0(2^) 
is a solution of the Dirac equation for a free particle. Since we assume that F'^'^ and are constants, the evolution 
operator can be expressed as a matrix exponential, 

C/(t,to) =exp|^7V^,(^^g'^(T-To)|, (3.9) 

that depends on constant four- vector q^^ satisfying the condition q^ — m^. 

Using ansatz p.8p we should choose such a basis in the space of solutions of the Dirac equation for the free particle 
that the action of the evolution operator on each element of this basis be reduced to multiplication by one and the 
same matrix function depending on quantum numbers of the basis. It is obvious that in our case we must choose the 
plane waves 

M^) = e^'^'^^l - Col'l^S^{q)){rq^. + m)V'o. (3.10) 

Here SQ{q) determines the direction of polarization of the particle; Co = ±1 is the spin projection on 6*0 (9); fpo is a 
constant four-component spinor. The wave function is normalized by the condition 



= m/qo, (3.11) 



and four-vector q^ is a kinetic momentum of the particle. Though the explicit form of a kinetic momentum oper- 
ator for a particle with spin interacting with dense matter and electromagnetic field is not known beforehand, the 
correspondence principle allows us to construct solutions characterized by its eigenvalues. 

Naturally, as far as equation (II. 3p is invariant under translations the canonical momentum operator = id^ is an 
integral of motion for this equation, too. However, the commonly adopted choice of eigenvalues of this operator as 
quantum numbers is not satisfactory if we prefer spin projection operators with clear physical meaning. The directions 
of canonical and kinetic momenta are different in the general case (see [23] ) and, as it was already mentioned in Section 
im projection of the spin is well defined in the rest frame of the particle where its kinetic momentum is equal to zero: 
q = 0. That is why in the construction of spin projection operators (see (12. 7p ) it is necessary to select unit vectors 
orthogonal to four- vector g'^. 
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Let us find the proper time t{x) and tlie phase factor K{x). Substitution of (|3.8p into (|1.3[) gives 



^^ioF'"'<y^.u - m| e-*^(^){/(r(x))*'o(a;) = 0, 



where = d^T. Since the matrix U{t{x)) is nondegenerate and the commutator [7^5^, ?7(r(a::))] is zero, we get 

rd,K{x) - ^7^/^(1 + 7') + ^7'f^M.¥''^9''7"^a - ^A^oi^'-V^, = 0. 

To solve this equation for iV^ and K{x), we should set the coefhcients at the linearly independent elements of the 
algebra of the Dirac matrices equal to zero. Thus 

d^K{x) = (3.12) 

and the system that defines the vector is 

VP^(m - (Nq)) + {Nip)q^ = 0, ^ioF^'^qa = -e^""^ N^^^q^, (3.13) 
The system is consistent provided that 

g^F'^V, = 0. (3.14) 
Equation (|3.14p must be held for an arbitrary q^, so we have 

*F^''F^, = ^^S^^F'^^F^^ = 0, (3.15) 

F^'•'f^ = 0. (3.16) 

Any antisymmetric tensor has an eigenvector corresponding to zero eigenvalue if and only if its second invariant 
I2 ~ \*F^'^F^^ is equal to zero. That is why conditions p. 151) and p. 161) are not independent and (|3.15p is a 
consequence of p. 161) . Comparing it with (|1.9p . we see that the condition obtained basing on the physical reasons 
alone and condition (I3.16P totally coincide. This ensures quasiclassical behavior of the particle in the background 
medium and allows obtaining the solution of equation (|1.3p in the form of ([33 
Using an orthogonal basis in the Minkowski space, 



1^ - /m 77^ - - — 77^ - — 77^ - ""^ ^ f3 17^ 

/A/ VA/" mVAfAf 

where Af = q^Hf^PHp^q" , Af = nlq^F^-P Fp^^q" , and relations (|B3]), (|B4| from AppendixEl we find 

= -q^^ ^liM + + 'll^ . (3.18) 

2{{Lpq)'^ — ui^tf^) 2(c/jg) 2{ipq){{Lpq)'^ ~ rrflLp^) 

According to the fact that f^^ and N^^ are constant values we obtain the proper time 

T = {Nx), (3.19) 

and the phase factor 

K{x) = {fx)/2, (3.20) 

which determines an energy shift of the neutrino in matter. 
Now we can derive the expression for the wave function: 

■^^Col^) = ^ E e-'^^^^Hl - Cl'l^S^M){l - Col'l.S^iqm'q, + H^o, (3.21) 
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where 



(3.22) 



SUq)^t^0^^£^. (3.23) 



System p.2ip represents the complete system of solutions of equation (II. 3p characterized by kinetic momentum of 
the particle q^ and the quantum number (^o = il which can be interpreted as the neutrino spin projection on the 
direction SQ{q) at r = {Nx) = 0. 

System p.2ip is nonstationary in the general case. The solutions are stationary only when the initial polarization 
vector (g) is equal to the vector of the total polarization 5'(|,(g) [ssj, S^{q) = S^^{q). In this case the wave functions 
are eigenfunctions of the spin projection operator ©tp — —l^l^iS^piji) with eigenvalues C, — ±1, and of the canonical 
momentum operator = id^ with eigenvalues P^. The orthonormal system of the stationary solutions, the basis of 
solutions of equation (|1.3p , can be written as (see Appendix |D| 

= e-(^c-)y'|J^(g)|(l - Cl'7^Si^^{q)){rq^ + m)^o, (3.24) 

where (q) is the transition Jacobian between the variables q^ and 

'dPl dPl dq°' 



Jf(g) = det(My) = det 



(3.25) 



dq^ dq^ dqi 

With the help of relations from Appendix |Bl the explicit form of the matrix Mij may be written as 

'^\m))v' 4(^g)((</^g)2-mV2)3/2 ' 

One can easily derive the following equality for arbitrary vectors g and h 

det ((5y + gihj) = 1 + (g • h). 

So we have 

'^^ ^ ^ 2yJ [fqY ^ m'^ip'^ ^ ^ ^ [fqY ^ TTi^f'^ ) 

Here Ii = jF^'^F^i, is the first invariant of the tensor P'"^. Note that to obtain a complete system of solutions for 
the antineutrino, it is necessary to change the sign of the kinetic momentum q^ . 

The dispersion law for the neutrino in dense magnetized matter is different from the one for the free particle and 
can be written as (see relations (|B7p . (|B8p and (|B9[) in Appendix [B|) 



P^ = rn? - /V4 - 2iilh - 2CA^ (Pl>)2 - $2^^ (3.27) 

where 

= P^ - /^/2, l>'^ = /''/2 + H'^'Pjm, 

A = sign (l + ^/.g^V/(2H-2M^/A ^ / ^ , UH^^h/m-A,ll, \ (3-28) 
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The appearance of the factor A in equation p.27p is a consequence of the fact that ( is projection of the particle spin 
on the direction defined by the kinetic momentum instead of the canonical one. 

In spite of the modification of the dispersion law described above, we see that the neutrino moving through dense 
matter and electromagnetic field may still behave as a free particle, i.e. its group velocity 

dP9 q 

V. = 4 = ^ (3-29) 

is the same for both polarization states of the particle. However, in interactions with other particles some channels 
of reactions which are closed for a free neutrino can be opened due to the modification of the dispersion law (see, for 
example, [H, H, 

Let us discuss now properties of nonstationary solutions in more detail. Solution (j3.21l) is a plane- wave solution of 
equation (|1.3|) . describing a pure quantum- mechanical state of a neutral particle with a nonconserved spin projection 
on the fixed space axis. Solutions (|3.2ip do not form an orthogonal basis. However, the considered system is not 
overcomplete, since the spectrum of the spin projection operator is finite. So the system can be easily orthogonalized. 
Generalization of the basis p.24p is 

= ^ E e-(^c-)y'|j^(^)|(i _ Cj'l,Si^M){l - Col\.S^{q)){rq, + m)^,,. (3.30) 
C=±i 

Thus we have just established that the unitary intertwining operator ()2.10p in our case is the Fourier integral 

operator [38] and it acts on elements of the plane- wave basis of solutions of the free particle Dirac equation p.lO|) in 
the following way: 

#,Co(^) = U{x,x,)M^) = i E e-*«^^-^)^)y^Jc(S(l - Cl'l,Si:^{q))Mx). (3.31) 

C=±i 

Action of the inverse operator is defined by the formula 

M^) = U-\x,xo)^,M - i E ^^"''^"'^'^^Tf^(l - C7Sp^rp(g))<^,Co(^)- (3.32) 

Since the intertwining operator is defined on the elements of the basis, its action on an arbitrary solution is defined 
as well. Hence, the explicit form of this operator as a function of coordinates and differential operators can be easily 
obtained. 

Note that an attempt to construct intertwining operator for equation (jl.3l) was undertaken in (39j for the case where 
only parameter is nontrivial. The result of the action of the operator suggested in ^ on plane- wave solutions of 
the Dirac equation for a free particle coincides with solutions obtained in our previous work for this particular 
case, and which can be derived from (|3.2ip . if one sets F^^^ — 0. We should emphasize that solutions p.2ip do not 
form an orthogonal basis. As a consequence of this fact, the intertwining operator suggested in '39'] is not unitary 
with regard to the standard scalar product (|2.6I) . 

The case of a massless neutrino is quite special. Equation p.7p obviously does not hold in the limit to -> because 
massless particle helicity is conserved, while the BMT equation describes spin rotation. To find wave functions 
of massless neutrino, one must take into account that from the mathematical point of view, the small group of 
representation of the Poincare group for a massive particle differ from that for a massless particle (see, for example, 
[25|). The small group for a massive particle is the rotation group of three-dimensional Euclidian space 0+(3), and 
nontrivial magnetic and electric moments are allowed. In contrast to that, the small group for a massless particle 
is the movement group of Euclidian plane E(2), which outlaws nonzero magnetic moment as well as electric one or 
other nontrivial coefRcients /i„ and e'„ in equation (jA2p . Therefore, one has to put /io = in equation (|1.3p . Then 
instead of (|3.24p . we obtain 

^,d^)U-o = e-'(«-)(i+^'')e-'(/-)(i-«)/2|i + ^,^1(1 „ C^5)^'^g^Vo, (3.33) 

for the neutrino, and 

'PMU^o = e'(«-)(i+f')e"'(/-)(i+0/2|i + ^^1(1 + C7')7^A.^o, (3.34) 
for the antineutrino. In these equations rj = 2{fq) ^'^'^ 1^ ~ ^- 
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IV. OPERATORS OF OBSERVABLES 



Let us find the explicit forms of the kinetic momentum operators and the spin projection operator ©tp in the 
coordinate representation. For this purpose we might exploit formulas (j2.12p : however, we follow a simpler way. 
The obtained solutions p.24p are classified by eigenvalues of the operators p'^ and Stp, so 



(4.1) 



Since solutions p.24p are also eigenfunctions of the canonical momentum operator p^^ = id^ with eigenvalues P^, we 
have 



(4.2) 



Now we should express eigenvalues of the kinetic momentum operator q^^ in terms of eigenvalues of the canonical 
momentum operator . From p.22p we have 



2{^q) 



i + C- 



(4.3) 



With the help of relations from Appendix |B] this can be rewritten as 



P2 _ ^2 + /2 /4 + 2f4li - ($/) 



The vector of total polarization in terms of the new variable is 



Si^Jq) 



($P)2 - m2|>2 



(4.4) 



(4.5) 



Here P, $ and A are given by (|3.28p . 



Because of (|4.2I) we can interpret P^f as a result of action of operator p^ 



id^ on the wave function. So by changing 



pfi ^ pfj. _ jiJ. j2 and /^/2 + i?'"'(py - /j./2) jm in formulas (|4.4p . (|4.5p . we obtain kinetic momentum operator 

and spin projection operator ©tp — ~7^1fj.Stp{q) in the explicit form. These operators are pseudodifFerential ones 
[38| and are determined on the solutions of equation (jl.3p with fixed mass m. 

To extend the domain of the definition of constructed operators, we need to replace mass m in (|4.4p and (14.51) by 
the matrix operator from equation (jl.3p : 



= ip- f/2f - f/4 + 2/^2/1 + 7V^VmP. + H^'l.IA^ + l') + 27'i/^''7MP- 



(4.6) 



In this way we obtain the covariant form for p*^ and ©tp. 

Unfortunately, the result of this substitution cannot be written as a compact formula, so we do not present it here. 
However, even if we do not know a covariant form of the operator p^, we may conclude that on the solutions of 
equation (|1.3p the relations, 



p2=TO2^ ^Mp^^m, (4.7) 

should hold. The first equality in (|4.7p is obvious; to proof the second one, recall that according to (j2.12l) 

TO = t/(x,xo)7>Mf^-'(a:,xo) = U{x,xo)rU-\x,x^)p^ - (7^ + P''((Z))p^. (4.8) 

From formulas (I3.3ip and (|3.32p it follows that r^{q) ^ S'(^(g), and as a consequence r^^{q)p^ — 0. So equation (|1.3p 
may be represented as 



(7'^p^ - m)W{x) 



(4.9) 
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Consider now special cases where the presented technique looks quite clear. Discuss the influence on the neutrino 
dynamics of the electromagnetic field alone, i.e. assume that = 0. In this case eigenvalues of the kinetic momentum 
operator would be expressible in terms of eigenvalues of the canonical momentum operator in the following 
way: 



Then the covariant form of the kinetic momentum operator is 



P^=P^ + ^==:^=6tp, (4.11) 



and the spin projection operator ©tp is defined by the formula 

Si-„ = siffn I 1 + 

Here 



6tp = sign ( 1 + ^ 6tp 1 6tp. (4.12) 



6.P ^ 7>^^^^- . (4.13) 

Thus operator Stp has a simple physical meaning. It characterizes a particle spin projection on the direction of the 
magnetic field in the rest frame of the particle. 

Vice versa, when the electromagnetic field is absent, but is nontrivial, we have 

and the covariant forms of the kinetic momentum and spin projection operators are 



ffj. ^fj. f2 _ fP-f f„] 



e.p = . (4.16) 

Note that if the matter is at rest and nonpolarized (f = 0), then 

6tp = sign(/")^^_P), (4.17) 



in other words ©tp is equal to the standard helicity operator up to the sign. 

We can find now spin projection operators for nonstationary wave functions p.2ip and p.30p . For this purpose, 
we introduce operators &± that act on the elements of system p. 241) as follows: 

<S+^,d^)^^^^^,i-0{x), &^W,^(x)^^^^^^^_^^{x). (4.18) 

Then operators ©i = ^(©+ + ©-), ©2 = 5^(6+ — ©-) and ©3 = ^©tp correspond to elements of the Lie algebra of 
the SU{2) group. Commutation relations for these operators are 

[©»,©j] =^ey■fc©fc. (4.19) 

To determine the explicit realization of operators ©± on eigenfunctions of operator let us choose the basis S'j^{q) 
(see (|2.8p ) in the form S^p{q), {q), S2 {q). Here S[^p{q) is defined by relation p. 231) . the other spacelike unit vectors 
are 

S>^( , ^o(g) + ^rp(g)('go(g)^tp(g)) ^ e^-P\.gop(g)5tpA(g) 



12 



As a result we have 



6± = 



1 V\Jc=±i{q)\ 
VI-^C=Ti(9)l 



where 



{Nx)\/ {(pq)"^ — ip^m? /m. 



(4.21) 



(4.22) 



Operators ©tp and 6± are integrals of motion. So the spin projection operator ©o that has eigenfunctions p.30p 
and eigenvalues Co = il is a linear combination of these operators. We can rewrite wave functions p.30p in the form 



So it is obvious that 



l-Co(5o(g)5tp(g)) 



n/2 



+ CoVl + Co(^o((?)5tp(9)) 



V2 ■ 



6o = -(5o(g)5tp(g))6tp + Jl - {So{q)S,p{q)f [©+ + 



(4.23) 



(4.24) 



Similarly to (|4.24p . one can construct the integral of motion ©o with eigenfunctions (|3.2ip and eigenvalues Co = ±1- 
Since 



^,Co{^) = Ji~C^{S,{q)SM) ^'i-'^''] ,1 + Cov'i + Co(^o(g)^tp(g)) ,'!f-'^^f 



we have 



©0 = -(5o(g)5tp((z))©tp + ^l~{So{q)SMY 
Or, in the other form, 



v/|Jc=-i(g)l 



©^ 



v/|^c=+i(g)l 
^/\Jc=-i{q)\ 



©_ 



©0 = 7'7p I (5o(g)5tp(g))^rp('z) - ^1 - (5o(g)5tp(q))2 [cos20 5r(g) - sin 20 5^' (g)] 



(4.25) 



(4.26) 



(4.27) 



Note that operator (|4.26p is not a self-adjoint operator with respect to the standard scalar product p.6p . It seems 
quite natural, since wave functions p.2ip do not form an orthogonal system. However, the system of wave functions 
is orthonormalized to the condition "one particle in the unit volume." In this sense wave functions (I3.2ip minimize 
the uncertainty relation which is due to (I4.19p 



((©l-(©l))2)((©2-(©2))^)--(©3)'. 



(4.28) 



Therefore, these wave functions describe spin-coherent states of the neutrino (on properties of coherent states see, for 
example, [401 )■ The given system of spin-coherent states is parametrized by four- vector 5'q (g). 



V. QUASICLASSICAL INTERPRETATION 



Let us discuss in more detail the physical meaning of spin-coherent states of the neutrino. For this purpose we 
construct vector and axial currents with the help of p.2ip . The vector current is 

= §{x)r^{x) - (5.1) 

and we see that (|3.2ip describes neutrino propagation with the constant velocity Vg^ — (l/q^ ■ The axial current is 

f^^§ixh^j>^^ix)^Co^S^\ (5.2) 

Here 

5'^ = -S^^{q){So{q)SM) + [S^{q) + S^^{q){So{q)Stp{q))] cos20 - -e'^''''\,Sop{q)Stpx{q) sm20, (5.3) 
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where 9 is determined by (|4.22l) . As expected (see p.Sp and p.6p ). vector and axial currents coincide with the 
solutions of the BMT equation, if the proper time is defined as in p.l9p . 

The three- vector of spin ^ can be expressed in terms of the four- vector components as 



C = s- 



(5.4) 



Then for C, we have 

C = Ctp (Co ■ Ctp) + [Co - Ctp (Co ■ Ctp)] COS 26 - [Ctp X Co] sin 29. 

Expression ()5.5p has a simple quasiclassical interpretation. 

The antisymmetric tensor G^'^ (see equation (|3.4p ) can be written in the standard form 



where 



G^'^ = (/xoP,/ioM), 



M = (/"q - q°f )/(2Mom), P = - [q x f ] /(2^o™) 



(5.5) 



(5.6) 



(5.7) 



Vectors P and M are analogous to the polarization and the magnetization vectors of medium. Note that the substi- 
tution F^^'^ F^^^ + G^'^ / jio implies that the magnetic H and electric D fields are shifted by the vectors M and P, 
respectively (we use here notation common in electrodynamics of continuous media [4l[) 



H^B = H-HM, D^E = D- 
Thus restriction ()1.9p in the explicit form is 

(E • f ) = 0, E/" - [B X f] = 0. 



(5.8) 



(5.9) 



This means that the Lorentz force and moment of force acting on matter are equal to zero, i.e. matter is at equilibrium 
state within the accuracy of our consideration. In particular, the vector of polarization of matter is parallel to the 
magnetic field if the matter is at rest. 

In the rest frame of the particle, equation p.3p can be written in the form 

C = 2^o[CxBo], (5.10) 

where the spin vector C is related to four-vector S^^ by equation (|5.4p and the value Bq is the effective magnetic field 
in the neutrino rest frame. This field can be expressed in terms of quantities determined in the laboratory frame, 



1 



>o = 
1 

TO 



,«B-[qxE]-4^-^) 



q"H - [q X D] - 



q" -f- TO 

q(q ■ H) ^ _q_ /^o 

q'^ + m 2/xo V 



(q-f) 

-f TO 



f 

2^0 ' 



(5.11) 



We see that the neutrino spin precesses around the direction Bo with the frequency lo — 2m|Bo|/g", the angle between 
Bq and the vector of spin being § — arccos(Co ■ Ctp)- The spin vector direction corresponding to stationary states Ctp 



is connected to the effective magnetic field as follows: 



Ctp — 



Bo 
IBnl 



(5.12) 



This fact explains in a simple way the stationarity of states with S^{q) = S'^^{q). 

Let us introduce a flight length L and a spin oscillation length Lqsc of the particle, remember that these are related 
hy 9 = ttL/Losc and that the scalar product {Nx) = t may be interpreted as the proper time of a particle. The spin 
oscillation length is 



27r|q| 



V(/g)2 - /2to2 - 4mH^-f^q, + AH^^'-H^^q^q- ' 



(5.13) 
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We can now write the probability for the neutrino which was arisen with polarization to change the polarization 
to —Co after traveling some distance L. This is 

= sin2?9sin2(7rL/Losc)- (5.14) 

It is clear that when (^g • ^^p) — Q ov 'd — 7r/2, the probability can be equal to unity and the resonance takes place. 
If the neutrino has the fixed helicity in the initial state 

S^{q) = ^{H,q\/H], Co = A, (5.15) 
m |q| 

then sin^ i? = 1 — (Bp • q)^/|BoP|qp and formula (|5.14p simplifies to the result widely discussed in the literature. This 
fact is not surprising. As was mentioned in Introduction, the Schrodinger-type equation with an effective Hamiltonian, 
which was used in papers 0] for obtaining this result, is merely the BMT equation in the spinor representation. 

VI. DISCUSSION AND CONCLUSIONS 

Let us compare our results with those of the standard quantum-mechanical approach to this problem. Since the 
neutrino behavior in dense matter under the infiuence of magnetic field, to the best of our knowledge, was not 
investigated before, we consider the neutrino propagation in a constant homogeneous magnetic field alone. 

In the studies of the influence of a stationary pure magnetic field on the neutrino spin rotation in the pioneer paper 
[isj , as well as in other papers, the stationary solutions 'Ppc,{x) first found in were used as the wave functions of 
a particle. These solutions are the eigenfunctions of the canonical momentum operator and of the spin projection 
operator Stp (see (14.121) ). The description of the neutrino spin rotation there is based on solving the Cauchy problem 
where the initial condition is chosen in such a way that the mean value of neutrino helicity is equal to ±1. It was taken 
for granted that the solution of the Cauchy problem can be expressed as a linear combination of the above-mentioned 
wave functions: 

w{x) = <^dp)i'vd^)- (6-1) 

C=±i 

However, such an assumption is incorrect. The point is that, once in a pure state the mean value of some spin operator 
is equal to ±1, then this state is described by an eigenfunction of this operator. In the general case, the construction 
of the eigenfunction of the spin projection operator as a superposition of only positive-energy solutions of equation 
(II. 3p is possible only when this spin projection operator commutes with the operator of the sign of the energy. The 
standard helicity operator (S • p)/|p| does not feature it. 

The given phenomenon is a sort of the famous Klein paradox [43| . To avoid the indicated difficulties, in relativistic 
quantum mechanics only self-adjoint operators in the subspace of wave functions with a fixed energy sign can be 
treated as operators of observables. The choice of integrals of motion as operators of observables is the necessary 
condition to satisfy this requirement [24]. 

In the case considered the canonical momentum operator is an integral of motion. However, the conserved operator 
of the spin projection which should set initial conditions to the Cauchy problem is uniquely — up to the sign — 
determined by the form of the Dirac-Pauli equation. This operator is ©tp- Therefore, it is impossible to construct 
a wave function describing a neutrino with rotating spin in the form of an eigenfunction of the canonical momentum 
operator for its arbitrary eigenvalues. The solutions similar to (|6.1I) can exist when the special values of the canonical 
momentum are chosen. So, if a neutrino moves parallel or perpendicular to a constant homogeneous magnetic field, 
eigenfunctions of the helicity operator are the superpositions of positive-energy solutions alone (4^ . 

To solve the problem we abandon the view that eigenvalues of the canonical momentum operator always impose 
a direction of the particle propagation. We found a self-adjoint operator which can be interpreted as kinetic 
momentum operator of the particle and obtained the complete orthonormal system of the solutions of equation (|1.3p 
with elements which are eigenfunctions of the given operator. On the base of this system we constructed solutions 
describing the neutrino with rotating spin. So our results enable one to treat a possible effect of the neutrino 
polarization change as a real precession of the particle spin. 

Consequently, the problem of neutrino spin rotation in dense matter and in strong electromagnetic fields is solved in 
full agreement with the basic principles of quantum mechanics. Using the wave functions of orthonormal basis p.24p 
or spin-coherent wave functions p.2ip . it is possible to calculate probabilities of various processes with the neutrino 
in the framework of the Furry picture. When choosing one or another type of the basis, it is necessary to take into 
account that, due to the time-energy uncertainty, stationary states of the neutrino can be generated only when the 
linear size of the area occupied by the electromagnetic field and the matter is comparable in the order of magnitude 
with the formation length of the process — the spin oscillation length in our case. 
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Appendix A 

Let us find the equation for the description of the neutrino behavior in dense matter and in electromagnetic fields 
at low-energy limit. When weak interaction with background fermions is considered to be coherent, the behavior of 
mass states of any one-half spin lepton should be described by the Dirac-type equation, 

(Z7^a^ + QJeni + 93„,attcr " m) W{x) = 0. (Al) 

In this equation the term 2Jcm describes interaction of the particle with the electromagnetic field and the term 2Jmattcr 
is responsible for weak interaction with matter. 

Following paper (45| . the nature of the interaction terms is determined by the restrictions that the equation be 
Lorentz covariant and gauge invariant; that the terms are linear in the electromagnetic fields and integral charac- 
teristics of matter, i.e. currents and polarizations of background particles; that terms do not vanish in the limit of 
vanishing momentum of the particle; that the charge and current distribution associated with the particle be suffi- 
ciently localized that its interaction with slowly varying electromagnetic fields and characteristics of matter may be 
expressed in terms of the electromagnetic and matter potentials and arbitrary high derivatives of these potentials 
evaluated at the position of the particle. These assumptions lead to the term Q3cm in the form 



n=0 



(A2) 



Here A^^ is the potential of the external electromagnetic field, F^'' — d^A" — d^A^^ is the electromagnetic field tensor, 
□ = d^dfi is the d'Alembert operator. The constants /ij^, e„, characterize the interaction, Eq is charge of the 
particle, fJ-o,£o are, respectively, anomalous magnetic and electric moments, and /i'q is an anapole moment. The 
expression for the QJmattcr can be found if we replace Afj_ in (jA2p with a linear combination of the currents (jl.l[) and 
of the polarizations (jl.2p of background fermions / with the proper choice of coupling constants. 

We have the minimal nontrivial generalization of the Dirac-Pauli equation neglecting terms with derivatives higher 
than the second in (jA2[) . In the expression for SJmatter we must hold only leading terms due to the proportionality of 
DA'^ to the sum of charged particle currents. Further restrictions for (|A1[) depend on the sort of lepton and on the 
model of interaction. The neutrino is a neutral particle, thus Eq = 0. In the framework of the standard model where 
it is assumed that the theory is T invariant and the neutrino interacts with leptons and quarks through left currents, 
its anomalous electric moment goes to zero (e'g = 0) and the term describing direct interaction of the neutrino with 
the currents contains multiplier (1-1-7^). As a result we come to equation (|1.3p with an effective four-potential 
which is determined by (|1.4p and (|1.5p . 

Appendix B 

For an arbitrary antisymmetric tensor A'"^, its dual tensor M^"^ = —^e'^'^P'^Apx, and for any four- vectors g^,h^, 
such as {gh) ^ the following relation takes place [46j : 

A''" {gh) = - [g'^A'^Php - A^'Phpg"] + * [h''*A''Pgp - ^A^Pgph"] . (Bl) 

This leads to the formula 

gp*A^^p*AP^g-h^ + {gp*A^^,h-f ~ hpAPpAP^h^g^ ~ [9^.AP,h-f = 
= gp-APp*AP^h-{gh) - g^A^ pAP ^h- [gh) = 2{ghfh. 

Here h = \Ai"' Ap^ = - \*AP''*Api, is the first invariant of the tensor A^" . 



(B2) 
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Appendix C 



The Lorentz equation for the four- velocity and the BMT equation for the spin vector 5''^ are 



m 

,2 



m 



(B4) 



Let Af"" = *Ft''', g^" = ip^", and h^' = q^' . Since *F^'•'Lp^ = -F'^'^ip^ = (see (pT^ . (IXTel) ). then from (HH) and 
(|B2l) we get 

2i/'^"^(^q) = if^^gp/'' - /^//""gp = miip^^r - /V'), (B3) 

and 

Let A^"" = 5^ = and h^' = P^'. Since ^F^''^^ = -F'^'^*^ = 0, then from dHT]), jHU we get 

2H^"'{^P) = H^^PPpf - f^'H^Ppp = m{^T - f^^l, (B5) 

and 

- = /2Ppff^"i?„.P'' + iUH^^^P.)^ = {.fP)UH''''H^uP'' - 2(/P)V^/i. (B6) 

It is possible to estabUsh by direct calculations using (jB4l) that 

(Pi') ^ (/!>) ^ / UH'^-q,/2m-2filh \ 
[qv) Up) \ ^{pqf - mV^ ; ■ 

It follows from ([B3| and (|B5l) that 

((^g) ~ ($p) ' 

and from jBl), dill), and dBTl) that 



(B8) 



(^<zr-mV = ^[(*Pf -m^l.^], 

($P)2 

(^9) ($P) (^9) 



(CI) 



Since u — 1 and {Su) — 0, we may rewrite these equations as [28| 

u'^ = f^'^'^u,., S"^ = O^^S-^, (C2) 

where 

r!^" = -P^" + 2^0(5"" - u''u")Fo,0{g^'' - u^u") (C3) 



is an antisymmetric tensor. 

From (jC2p it is obvious that the evolution operators for the Lorentz and BMT equations are the same and fulfill 
the relation 

P^''(r,ro) = f^^„P"''(r,ro). (C4) 
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From (|331) and (|3^ it follows that 

j(;(t) = i?^,(r,To)jj:.(To) = .^(t)7^*'(t) = l^(To)f/(To,T)7^C/(T,To)tf'(To), (C5) 

and 

Therefore, the equation for the evolution operator of the quasiclassical spin wave function takes the form 

?7(r,ro) -Z{/(r,ro), (C7) 

where Z must obey the relations 

[7^Z] = ^]'^^7., [7V,^] = f^^V7.- (C8) 



It is obvious that 



Using the relations 



and 



Z = -n^.a'^r (C9) 



^F^'^^M- = «7'*^''"^.7p^"7a + F'^'^u^^^^u'^^^, (CIO) 



we find 



F^V^, = z7'*i^''V^,, (Cll) 

Z = i7' (£^^'''7^.7. + MoJ^^"M,7^7""a) • (C12) 

We must take into account that the electric charge of neutrino e is equal to zero. Replacing ^^^F^" to fio*F'^^ + 
{f^u'^ — u^f^)/2 (see l|3.3p and p.4p ). and introducing the notation = m,u^ we obtain equation (|3.7p . 

Appendix D 

Let us prove the system (I3.24|) is orthonormal, i.e. 

= j d^^l,^,{x)1',c{x) = {2T:f5\ci-c()5a' (Dl) 
for the solutions with the same signs of energy and 

= j d^^l,^,{x)^q^{x) = (D2) 
for the solutions with the different signs of energy. After integration we get 



= (27r)^yGc^/7c(±^(^±(<z')) (g)'5'(Pc'(±9') " ^ciM)), (D3) 
N^. = {2^)\lj^.{±q')Jj^{Tq) (^±(g'))%r(9)'5'(Pc'(±9') ^ Pc(Tg)), (D4) 



where the conventional eigenspinors tp^{q) obey the equations {'^^q^ T = 0; J^((j') is the Jacobian for 

transition between the variables q^^ and P^{q). 
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In the standard representation for the gamma matrices, we obtain exphcitly [47 



y/2qf^{q" +m) V ('^ ' 

{a- ■ q) u}7 



where ct; are the PauH matrices and factor 1 / A/2g0(g° + m) has been included for normahzation p. lip . These spinors 
satisfy the orthogonahty relations 



[i;f{q)) ^ 7°# (g) = 0, [^f (g)) ^ ^f, (q) = ^cC , 



(D5) 



if Lu^ are nonvanishing, but otherwise arbitrary two-component spinors which are chosen such that (w^)tw^ = '^CC- 
It is convenient to use the remaining uncertainty in co^ to require ipfiq) to be eigenstates of the spin projection 
operator 



-l'l,S^piq)i't{q)^C^bf{q). 
For this purpose ui^ should be eigenspinors of three-dimensional spin projection operator (see (|5.4I 



(D6) 



(D7) 



To calculate spinors lu^ we have an opportunity to choose a special reference frame. Let us take the reference 
frame where = {/, 0,0,0}, H = (0,0, -ff), and E = 0. It is possible, if is timelike four-vector and the relation 
Ft^'^f^ = is fulfilled. 

In this reference frame the vector of the total polarization is 



Sf^iq) = ^^^{^-L + qiiq±/q±, q°cos(f>, q°sm<j), q°q±/q±Y 



(D8) 



and 



q±R± 



It 



/(gO + m) 



± 2//oH ■ 



(D9) 



where 



i?± = sign(/)./(l + AfilHyp) + {q±lqj_f 



(DIO) 



Here we use the notation 



\/{q^Y + {q^Y = q±, ^ q± cos(f), q^ = q±_ sin0, q^ = g,,, q^^±2nomH/f = q±. 



(Dll) 



The explicit form of spinors uj^ is 



(gO + to) uj^ = - 



e-#/2 [(qO + m)T2fioHqJf] 
„»0/2 [(gO ^ ni){±CR± - q±/q^) T 2fioHq^/f] 



ia-Cl)u^^^^ [q^ q, (±Ci?± - q^/q^)] ) ' 



(D12) 



where 



T± = 2R± 



qi- f{q° + m) 



1/2 



(D13) 



19 

In the selected reference frame the components of canonical momentum for positive-energy and for negative-energy 
solutions, respectively, are 

p,(±g)-±gxcos0fi + C ^^^^^ 

P^i±,) = sin0 (l + C^^l;^ j , (D14) 

Consider at first the C' = C case. In this case on the right-hand side of equation (ID3p we have 

5'{Pd±q') - PdM)) = I Jc(±g)r^<53(q' _ ^15) 
and as far as q'^ = = we get q'^ = g^. Therefore, 

= (27r)353(q_q'). (D16) 

In equation (jD4[) the delta function provides the relation 

q + q' = 0. (D17) 

So we have 

= 0. (D18) 
Consider next the C' 7^ C case. The delta function provides the following relations: 

q'±/q'i_ = q±/q±, = (1^, (D19) 
on the right-hand side of equation (|D3|) and relations 

q'±/q'^ = - q^/q^, 0' = + ^, (D20) 



on the right-hand side of equation (jD4p . 

Using (jD19P and ()D20p and bearing in mind that (' = we find 



(V'^^(g'))%f('7) = 0, (V.±^(<z'))%c^(g)=0. (D21) 

Thus we have proved relations (jDl[) and (jD2[) . i.e. the orthogonality of system p.24p . It is easy to verify that for 
the cases where four- vector is spacelike or lightlike we can get the same result. 
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